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Abstract 

We show in the framework of a tractable model that revivals and fractional revivals of wave 
packets afford clear signatures of the extent of departure from coherence and from Poisson statistics 
of the matter wave field in a Bose-Einstein condensate, or of a suitably chosen initial state of the 
radiation field propagating in a Kerr-like medium. 
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The evolution of a quantum wave packet subject to a confining potential is of current inter- 
est in several experimentally realizable situations. Foremost among these is the dynamics of 
the matter wave field of a Bose-Einstein condensate (BEC) confined by a three-dimensional 
optical lattice The condensate is, in general, in a coherent superposition of different 
atom-number states ODD with a repulsive interaction between the atoms. This state 
evolves in time in the confining potential. If the number of atoms and the number of lattice 
sites are both sufficiently large, the atom number distribution in each well obeys Poisson 
statistics to a good approximation. With increasing well depth and decreasing tunneling 
energy, the wells can be taken to be sufficiently isolated from each other. The atom number 
distribution in each well then departs from the Poisson, and significant non-classical effects 
manifest themselves jl]]. These include squeezing and dissipation, sub-Poisson statistics, as 
well as revivals and fractional revivals at particular instants of time. A realistic and tractable 
model for the dynamics of the atoms in each well is provided by the Hamiltonian 

H = hxa ]2 a 2 = h X N(N - 1), (1) 

where a , a< and N are annihilation, creation and number operators of atoms, and \ charac- 
terizes the energy needed to overcome the inter-atomic repulsion in adding an atom to the 
population of the potential well. 

As is well known, the formal equivalence between bosonic atoms and photons has been 
exploited in recent years to bring out the deep analogies between quantum optics and atom 
optics involving BECsjsj], enabling a fruitful two-pronged approach to the problems of quan- 
tum computing. Wave packets propagating in a nonlinear optical medium display precisely 
the variety of non-classical features [6] mentioned above that BECs display. In particular, 
revivals and fractional revivals, which are now recognized to be generic features of wave 
packet evolution in nonlinear quantum dynamics, constitute^ striking aspect of experimen- 
tal observations [l o 
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I12I l^j j. including that of the dynamics of a single- mode field 



propagating in a Kerr-like medium. In this case the initial wave packet is a coherent state 
of the radiation field, and it is precisely the Hamiltonian in Eq. (pj (where the operators 
now refer to photons) that governs the dynamics. 

The identification of clear signatures of revivals and fractional revivals helps distinguish 
between wave packets that obey Poisson statistics and those that obey sub-Poissonian statis- 
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tics, and also provides valuable information on the degree of coherence enjoyed by the system. 
This would also be of practical importance in quantum computation using wave packets, 
where logic gate operations are envisaged to be implemented at the precise instants of frac- 
tional revivals 3] • 

Specifying the state of a BEC in an actual experiment is not simple, and several models are 
extant. However, plausibility arguments may be given to support a pure state description Q] 
of the BEC, according to which the state at any time t has the general form 

W0> = £^|n>, (2) 



where the expansion coefficients b n (t) are model-dependent. It is of great interest to study 
the departure from coherence of an initial state under time evolution governed by a hermitian 
but nonlinear Hamiltonian as in Eq. (0). What is required for this purpose is a model initial 
state that has three basic features: a precisely quantifiable, preferably tunable, degree of 
departure from perfect coherence, sub-Poissonian statistics (a standard deviation that is less 
than the mean), and phase-squeezing. All these properties are possessed by the normalized 
state 

(qt)"» (qt)"|q) 
\a,m) = — . = — . = (3) 

y/(a\ a m a) m \a) y/m\L m (-u) 

where m is a positive integer, \a) is the standard oscillator coherent state defined by a \a) = 

a \a) , a G C, v = \a\ 2 , and L m is the Laguerre polynomial of order m. The departure of 

\a,m) from perfect coherence arises due to the addition of m atoms to \a), a feature that 

becomes more pronounced with increasing m. In the context of quantum optics, in which 

this state has first been studied 15j], \a, m) is called an m-photon-added coherent state, and is 

produced in laser-atom interactions under appropriate conditions. Its non-classical features 

include both phase squeezing and sub-Poissonian statistics. The latter property implies that 

the standard deviation in the atom number N of the state behaves like N^~@ rather than 

iVa, the exponent f3 being a calculable decreasing function of m. 

While \a,m) is not an eigenstate of a, it may be regarded^] as a "nonlinear coherent 

state" , in the sense that it is an eigenstate of the operator [1 —m a with eigenvalue 

a. The state \a, m) can also be viewed in another way. Instead of the Fock basis {|^)}, we 

may consider the unitarily transformed basis {|n, a)} formed by the generalized coherent 

states |n, a) = exp (aa^ — a* a) \n) . ( Equivalent ly, for a given n, \n, a) is simply the state 
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(a* — a*) n \a), normalized to unity.) The initial state |^(0)) of the condensate can be 
expanded in the basis {\n, a)} instead of the Fock basis. Likewise, it can be shown that, for 
a given m, the state \a, m) is a superposition of the form Yl™=o c n \ n i a )- Thus, in practice, 
\a,m) is a very appropriate candidate for the initial state of the condensate. 

In this paper we show that distinctive signatures of revivals and fractional revivals of 
the condensate are manifested in the mean values of certain basic operators pertaining to 
the system, and in their variances. We examine the precise manner in which the departure 
from coherence of the initial condensate wave function affects its subsequent dynamics, 
particularly at the instants of revivals and fractional revivals. The distinctions between 
different fractional revivals that occur in between two successive revivals are also brought 
out. 

The matter wave field in a BEC is essentially given by the expectation value (i()(t)\a\ip(t)) . 
Its real and imaginary parts are the expectation values of the hermitian combinations (a + 
a')/2 and — i(a — a')/2, which in turn correspond to the cases ip = and <p = —\ts of the field 
quadrature £ = [ae lLp + a) e~ %v )/2 that is customarily taken 15] as the basic observable in 
the quantum optics context. We therefore set x = (a + a^)/2 1//2 and p = —i(a — a^)/2 1 ^ 2 (so 
that [x, p] —i), and examine the expectation values and variances of x and p as the system 
evolves from the initial state in Eq. Q under the Hamiltonian H of Eq. (£[]). We shall see 
that the time dependence of these mean values and the corresponding variances mirrors, in 
distinct ways, the occurrence of different fractional revivals between two successive revivals 
of the initial state. We use the convenient notation 

(x{t)) m = (a, m\e iHt/h xe- iHt/h \a, m) , (4) 

with an analogous definition for (p(t)) . As H is diagonal in the number operator N, it 
follows that the mean atom number, given bvjl5| 

/Art (m + l)L m+1 (-u) 

remains constant in time. When v = 0, (N) reduces to m, as required; while for v 3> 
m, (N) = v + 2m + 0(v~ l ). The moments of N, and hence the sub-Poissonian statistics 
of the atom number, also remain unaltered in time. 

It is helpful to use as a reference, for the purposes of subsequent comparison, the results 
that obtain in the case when m — 0, i.e., for an initial state that is just the coherent 



state | a) (which has, of course, a Poissonian number distribution with mean value v). A 
straightforward calculation gives 

(x(t)) = exp [-u(l -cos 2yt)\ 

x {x cos (z/sin2xt) +Po sin (z/sin2xt)}, (6) 

(p(*))o = ex P [ _z/ ( 1 _ cos2 X^)] 

x {po cos (z/sin 2%i) — xo sin (V sin 2;\i) } • (7) 

Here a = (x Q + ip )/2 1 ^ 2 so that v = \a\ 2 = \{x$ + pi) ■ The parameters x and p signify 
(in the case at hand, namely, for m — 0) the locations of the centers of the initial Gaussian 
wave packets in position and momentum space, respectively. It is evident that (x(t)) and 
(p{t)) are periodic in t, with a period it/x = T Tev , the revival time. (For the Hamiltonian 
in Eq. T rev coincides with the "classical orbit time" T c \ because the coefficients of the 
terms linear and quadratic in N happen to be equal in magnitude.) In fact, at the level of 
expectation values, this case can be cast in the form of a classical nonlinear oscillator with 
"dynamical variables" X = (x(t)) exp [i/(l — cos 2%t)] and Po = (p(t)) exp cos 2%t)] : 
setting Zo(t) = exp {iv sin 2%t), Eqs. (0) and (J2J) become 

Xo(t) = x Q Rez (t) +p lmz (t), 

Po{t) =p Rez (t) -x 1mzo(t). (8) 

If we now re-parametrize time according to r = sin 2%t, we have dX /dT = uP Q and 
dP^/dr = —uXq. Thus, formally, we have essentially a nonlinear oscillator of frequency 
|[X^(0) + P 2 (0)], X o (0) and Pq{0) being the initial values of X (t) and Pq{t), respectively. 

The time dependence of (x(t)) and (p{t)) for m^fl differs in striking ways from the 
foregoing, even for small values of m. The revival time T rev remains equal to it/ x, of course, 
but in the intervals between revivals the time evolution is considerably more involved than 
the expressions in Eqs. (JBJ) and (JZJ) for the case m = 0. This implies that even a small 
departure from coherence in the initial state of the condensate and from Poissonian number 
statistics leads to a very different time evolution of the system and the phase squeezing it 
exhibits. We have calculated the exact expressions for (x(t)) m and (p(t)) m , and these are 
as follows. Their initial values are given by 



where Lm(—v) = dL m+ i(— v)/du is an associated Laguerre polynomial. Analogous to 
(X , P ) , let us define 

X m (t) = (x(t)) m exp [u (1 - cos 2 X t)\, 

Pmit) = (p(t)) m exp [u (1 - cos2 X t)]. (10) 
The solutions for these quantities may then be written in the compact and suggestive form 

X m (t) = x Rez m (t) + p lmz m (t), 

Pm(t) = p Rez m (t) - x lmz m (t), (11) 

where 

L^(—ve 2ixt ) 

Zm(t) = \ — ex P [* ( 2m Xt + v sin 2 X t)} . (12) 

A number of differences between these results and those for the case m = are note- 
worthy. First, \z m \ varies with t, in contrast to |#o(£)l = 1- (Note also that z m (0) = 
Lm (—u) I L m (—v) 1.) The time dependence of X m and P m involves the sines and cosines 
of the set of arguments (2\ It + u sin 2\t) where I = m, . . . , 2m. Thus, not only are "higher 
harmonics" present, but the arguments also involve secular (linear) terms in t added to the 
original v sin 2\t. This important difference precludes the possibility of subsuming the time 
dependence into that of an effective nonlinear oscillator by means of a re-parametrization of 
the time, unlike the case m — 0. 

We present the rest of our results with the help of figures based on numerical computation. 
In what follows, we set x = 5 for defmiteness, and also restrict ourselves to the case Xq = po 
(there is no significant loss of generality as a result of this symmetric choice of parameters) . 
The presence of the overall factor exp \—v (1 — cos2xt)] in the expressions for {x(t)) m and 
(p(t)) m implies that, for sufficiently large values of the parameter u, the expectation values 
remain essentially static around the value zero, and burst into rapid variation only in the 
neighborhood of revivals. Smaller values of v enable us to resolve the details of the time 
variation more clearly. 

FiguresIHa) and (b) are, respectively, plots of the expectation values (x(t)) and {x(t)) 10 
versus t (in units of T rev ) for parameter values xq = po = 1 (i.e., for v = 1). The revivals at 
integer values of t/T rev are manifest. With increasing m (or a decreasing degree of coherence 
in the initial state), the relatively smooth behavior of (x(t)) gives way to increasingly rapid 
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FIG. 1: (x(t)) q as a function of time 
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FIG. 2: "Phase plot" of (p) us (a) 



oscillatory behavior in the vicinity of revivals. The range over which the expectation value 
varies also increases for larger values of m. Essentially the same sort of behavior is shown 
by (p(t)) m . However, a "phase plot" of (p(t)) m versus (x(t)) m in Figs. Efa) and (b) reveals 
complementary aspects of such oscillatory behavior with increasing m, showing how the 
oscillations in the two quantities go in and out of phase with each other. (The entire closed 
curve in each case is traversed in a time period T rev .) 

The initial state also undergoes fractional revivals in the interval between any two succes- 
sive revivals. In principle, it can be shown that, in the interval (0, T rev ), fractional revivals 
occur at the instants T rev /k where k = 2, 3, ... , and also, for any given k, at the instants 
j T rev / k where j = 1,2,... , (k — 1) . These fractional revivals are signaled by the appearance 
of k spatially-distributed wave packets similar to the wave packet representing the state at 
t = 0. The fractional revivals at the instants jT Tev /k show up in the rapid pulsed variation 
of the k th moments of x and p, and not in the lower moments[17]. However, if we use \a, m) 
as the initial state, then, even for relatively small values of m, the signatures of fractional 
revivals appear for values of v that are not large, in contrast to what happens when the 
initial state is the coherent state \a) . (Recall that (N) is determined by v according to 




FIG. 3: Variation of the uncertainty product with time 

Eq. ©.) 

For illustrative purposes we investigate the specific case of the fractional revival at t — 
|T rev . This corresponds to the appearance of two spatially separated similar wave packets, 
i.e., a single qubit in the language of logic gate operations. Plots of the product Ax Ap of the 
standard deviations of x and p versus t/T rev over a full cycle are shown in Fig. El for initial 
states given, respectively, by the coherent state \a) (dotted curve), the atom-added state 
\a, 1) (dashed curve), and the multi-atom-added state \a, 10) (bold curve), for the same 
parameter values as above (x = 5, x$ = po = 1, so that v = 1.) It is seen that hardly any 
trace of the fractional revival is evident in the case of an initially coherent state, in marked 
contrast to the case of the atom-added states, in which the fractional revival is signaled by 
oscillations whose frequency and amplitude increase quite rapidly with increasing m. This 
effect gets masked for larger values of the parameter z/, when these oscillations are relatively 
insensitive to the value of m. 

Another striking feature that provides a clear distinction between the revivals at t = n T rev 
and the fractional revivals at t = {n + ~) T rev is illustrated in Fig. 0] (a), which is a plot 
of Ap versus Ax . The dotted and full lines correspond to m = and m = 5, respectively. 
At t — 0, these quantities are equal, and have small values. As t increases, they rapidly 
build up, oscillating about the radial Ap = Ax line with an initially increasing, and then 
decreasing, amplitude. A maximum value of Ax and Ap is attained, at which these quantities 
then remain nearly static, till the onset of the fractional revival at T rev . They then begin 
to oscillate rapidly once again, but this time in a tangential direction, swinging back and 
forth along an arc with an amplitude that initially increases and then decreases to zero: 
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FIG. 4: "Phase plots" of higher moments of p and x 

in other words, the individual standard deviations fluctuate rapidly in the vicinity of the 
fractional revival (while [(Ax) 2 + (Ap) 2 ] 1 ^ 2 remains esentially unchanged in magnitude), in 
marked contrast to what happens at a revival. It is evident that all these features are very 
significantly enhanced and magnified for non-zero values of m, relative to what happens for 
the case m — 0. 

Similar signatures of fractional revivals for higher values of k can be discerned by using 
initial states \a,m) even with relatively small values of u, the oscillations in the moments 
of observables becoming more pronounced with increasing m. For instance, signatures of 
the fractional revivals at t = |jT rev , j = 1, 2, 3 are clearly discernible in the behavior 
of the fourth moments of x. Figure 0] (b) is a plot of ({Sp) 4 ) versus ((5x) A ) where 5x = 
x — (x) , 5p = p — (p) , for v — 1. The dotted and full lines correspond to the initial states 
|a) and |a, 2) respectively. Once again, the magnification of the variations that occurs for 
even a small value of m is manifest. 

In conclusion, we have demonstrated that an atom-added initial state of the form \a,m) 
shows significantly increased sensitivity to revivals and fractional revivals. In this sense, one 
may therefore expect that the inevitable departure, in practice, of the initial state of a BEC 
from perfect coherence can in fact be used to advantage. 
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